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Abstract We study theoretically the dynamics of charged quantized vortices (CQVs). 
CQVs (ion-vortex complexes) have been used as an important probe in the field of 
, , quantum hydrodynamics. Recent experimental studies of quantum turbulence at very 

low temperatures utilize CQVs. In this work we propose the equation of motion of 
CQVs based on the vortex filament model. An analytical solution for a charged vortex 
ring shows how it expands under an electric field. Numerical simulation reveals the 
characteristic behavior of CQVs under an electric field. 
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1 Introduction 

^ ■ CQVs (ion-vortex complexes) have been used as an important probe in the field of 

quantum hydrodynamics since the pioneering work by Rayfleld and Reif pQ. Recent 
experimental studies on quantum turbulence at very low temperatures are performed 
by means of CQVs [2]. However, there has been very little theoretical information on 
the dynamics of CQVs [3JH]. In order to understand the recent experiments, we need 
ON ' the information of the nonlinear dynamics of CQVs under an electric field. In this work, 

first, we propose the equation of motion of CQVs on the vortex filament model. Then 
■ we discuss some typical analytical and numerical solutions. 



2 Equation of motion of CQVs 

First, we remember how to derive the usual Schwarz's equation [5] 

^=v s + as'x (v n - v s ) - as x [s x (v n - v s )] . (1) 
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Here s refers to a point on a vortex filament and s' is its derivative with respect to 
the coordinate along the vortex. The symbols v s and v n are the superfluid and the 
normal fluid velocity fields, and a and a are the coefficients of the mutual friction. 
The Magnus force per unit length is written as 

fM = PsKS x - v s ) (2) 

with the superfluid density p s and the circulation quantum K. The drag force per unit 
length due to the mutual friction is 

f D = -ap s KS X [s X (v n - V s )] - a p s KS X (v n - V s ) . (3) 

Since the inertia of the vortex core is negligible, the equation of motion per unit length 
is fjvf + $d — 0, which leads to Eq. (pQ). 

Next, we consider the case that a part of a vortex is electrically charged. This 
describes an ion (the charge e and the radius R[ on ) trapped by a vortex core, namely 
an ion- vortex complex. An electric force f# = eE acts on the charged part. Now the 
equation of motion is 2R- lon (i^[ + fp) + f# = 0, from which we obtain 

"77 = v * + 7T75— s ' xE + as'x(v n -v s )- as x [s' x (v n - v s )] . (4) 

at ZK[ on p s K 

The second term of the r.h.s. is characteristic of the present case. This term is explicitly 
included in the vortex filament model for the first time. The physics of this term has 
been already discussed in some previous literature [1,3, 4 , which was limited to simple 
symmetric cases. In the rest of the paper we consider the case of zero temperature 
without mutual friction. 



3 Analytical solution of expansion of a charged vortex ring 

A key point of this issue is how the ion is localized along a hollow vortex core. Early 
works speculated that the ion is in effect distributed uniformly along a core PQ. The 
observations of the ion mobility along the vortex core show, however, that the mobility 
increases rapidly with decreasing temperature [6 j. The physics of the behavior is not 
yet fixed. In this work we consider both cases of the uniformly distributed ion and the 
localized ion, not confining ourselves to either model. 

The model of the uniformly distributed ion enables us to obtain the analytical 
solution of how a charged ring expands under an electric field. Now we assume that 
initially a charged vortex ring (CVR) with a radius Ro lie on a plane perpendicular to 
an electric field E in absence of any flow (v n = v s =0). The charge e is supposed to 
be uniformly distributed along the ring. The equation of motion J4} of a ring with a 
radius R is reduced to 

dt 2ttpskR 

whose solution is 
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Fig. 1 Configuration of a vortex loop at t = 0.15 msec. These three figures depict the same 
vortex from different directions. The electric field E is applied along the z axis. The small 
solid circle, located at rightmost extent of the vortex in each of the three frames, shows the 
localized ion. The arrows on the vortex lines show the direction of circulation. The z axis is 
exaggerated by a factor of 500. The ring radius in (c) is 0.5053 /xm. 



4 Numerical simulation of CVR dynamics 



We performed the numerical simulation of a CVR dynamics by the vortex filament 
model [7 , and found some behavior characteristic of the system. The simulation begins 
with a perfectly circular vortex ring. The ion is assumed to be localized on a point on 
the vortex ring. The ion is pulled by an electric field, exciting Kelvin- waves (KWs) with 
small amplitude. The KWs propagate along the vortex ring to keep approximately the 
circular ring shape. As a result, the expansion of the ring with a localized ion is almost 
same as that with a delocalized ion. 

Consider a CVR with radius R — 0.5 /xm whose self-induced velocity is towards 
the +z direction. An applied electric field is Ez with E — lOOV/cm, and the ion 
radius is Ri on = 15 nm. These parameters follow the recent experiments by Walmsley 
and Golov [8]. We made the numerical simulation by the vortex filament model. The 
motion of a charged part is described by Eq. J2}, while other neutral parts move with 
the usual equation neglecting the second electric term. The numerical space resolution 
is A£ — 3 nm, and the time resolution is At — 0.1 nsec. 

Figure 1 shows the vortex configuration at t = 0.15 msec. Since the deformation 
along the z axis is quite small, it is exaggerated by a factor of 500. Since the self- 
induced velocity of the ring is along the electric field, the charged part of the ion likes 
to expand through the second term of Eq. J4}. The resulting cusp excites KWs which 
propagate over the whole ring. However, these KWs are never amplified too much, so 
that the vortex ring propagates keeping the circular shape. This is consistent with the 
analysis based on the mode expansion |4]. 

Figure 2 shows the comparison of the radius growth between the solution of Eq. JB} 
for the delocalized-ion model and the numerical simulation assuming the localized-ion 
model. The agreement is quite nice. This is because the KWs starting from the ion 
soon cover the whole ring, as discussed in the next section. Although the observed time 
period is short, we should expect this agreement for longer time too. As the vortex ring 
becomes bigger, the propagation velocity is reduced. 
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Fig. 2 Expansion of the ring radius. The solid line shows the numerical result for the localized- 
ion model, while the dashed line refers to the analytical solution of Eq. ((6} for the delocalized- 
ion model. 



5 Discussions: the role of KWs 

It is well known that KWs play important roles in quantum turbulence at very low 
temperatures [9 , 10 . So do they for the present problems. The key question is how KWs 
transfer the energy from the work due to the electric field to the whole vortex ring to 
expand it. 

The early elementary process is clearly shown in Fig. 3. The electric field pulls 
the ion to excite a small cusp there. The resulting KWs have initially the wavelength 
comparable to the cusp scale, propagating apart from the ion. The nonlinear interaction 
is known to work between different wavenumber modes of KWs and cause the KW 
cascade [11]. As the initial KWs propagate apart from the cusp, KWs with smaller 
wavenumbers are also excited through the nonlinear interaction to follow the initial 
KW. Figure 4 shows how the KWs generated from the ion propagate and cover the 
whole ring. Through this process the energy injected from the ion becomes about 
uniformly distributed over the ring, so that the vortex keeps the ring shape. This is the 
reason why the radius growth of the solution of Eq. J6} for the delocalized-ion model 
agrees with that of the simulation using the localized- ion model (Fig. 2). 

One mystery is why the KWs are not so much amplified in spite of the continuous 
excitation by the electric field. One reason would be that the energy of KWs is partly 
consumed in expanding the vortex ring. Another reason may come from the finite 
numerical space resolution. The simulation cannot describe the KWs with wavelengths 
shorter than the resolution and this works in effect as some dissipative mechanism, 
which occurs in the simulation of the KW cascade [11] , 

In order to justify quantitatively the arguments in this section, we need the analysis 
of the spectra of KW during the dynamics, which is now in progress. 



6 Conclusions 



We study theoretically and numerically the dynamics of charged vortex rings under an 
electric field. First we introduce the characteristic term into the vortex filament model. 
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Fig. 3 Early elementary process of how the cusp grows from the ion and Kelvin-waves are 
excited. The ion is at the top of the cusp. The electric field is applied along the z direction. 
The time is 0.7nsec, 1.7nsec, 4.0nsec, 7.0nsec and lO.Onsec from the left to the right. The z 
axis is exaggerated by a factor of 500. 
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Fig. 4 Some late stage showing how the Kelvin- waves excited at the ion propagate and extend 
over the whole vortex ring. The solid circle refers to the ion. The time is 30nsec, lOOnsec, 
200nsec, 400nsec and lOOOnsec from the left to the right. The z axis is exaggerated by a factor 
of 500. 



Next the dynamics of a charged vortex ring is investigated by two models. One model 
assumes that the charge is uniformly distributed over the vortex because of the high 
mobility of the ion. Then we obtained the analytical solution of how the ring expands. 
The other model supposes that the ion is localized on the vortex. The vortex filament 
simulation reveals how the electric field causes a cusp at the ion and the KWs are 
excited to extend over the vortex. 

These findings would be helpful for understanding the recent experiments of quan- 
tum turbulence using ions by Walmsley and Golov [8 j. It should be quite interesting 
and important to study the dynamics of multi-charged vortices and charged vortex 
tangles as studied by them. These developments are now in progress. 
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